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We estimate the quantum corrections to the ground state 
energy in neutron matter (which could be termed as well ei- 
ther shell correction energy or Casimir energy) at subnuclear 
densities, where various types of inhomogeneities (bubbles, 
rods, plates) are energetically favorable. We show that the 
magnitude of these energy corrections are comparable to the 
energy differences between various types of inhomogeneous 
phases. We discuss the dependence of these corrections on 
a number of physical parameters (density, filling factor, tem- 
perature, lattice distortions). 

PACS numbers: 21.10.Dr, 21.65.+f, 97.60.Jd 

The investigation of the nuclear matter in the neu- 
tron star crust below the saturation density leads to the 
consideration of exotic shapes of the nuclei immersed in 
a neutron gas. It was realized long ago Q| that when 
the nuclei in dense matter occupy more than half of the 
space it is energetically favorable to "turn the nuclei in- 
side out" and form a bubble phase ■ To date a large 
number of calculations have been performed pertaining 
to the structure of the neutron star crust. The liquid drop 
model or the Thomas-Fermi approximation calculations 
|^| |llj predict rather small energy differences between dif- 
ferent phases, of the order of a few keV/ fm 3 . (N.B. even 
though we often refer to energy, we actually mean energy 
density.). Apparently an agreement has been reached 
concerning the existence of the following chain of phase 
changes as the density is increasing: nuclei — > rods — > 
plates — » tubes — > bubbles — > uniform matter. The den- 
sity range for these phase transitions is 0.04 — 0.1/to~ 3 
P,Eo|. Moreover, it was established that these phases 
exist up to temperatures of about lOMeV ||. At den- 
sities of the order of several nuclear densities the quark 
degrees of freedom get unlocked and the formation of var- 
ious quark matter droplets embedded in nuclear matter 
becomes then energetically favorable [Q. 

The appearance of different phases is attributed to 
the interplay between the Coulomb and surface ener- 
gies. Since most of the published works were based on 
the minimization of some density functional in a single 
Wigner-Seitz cell, the calculation of the shell correction 
or Casimir energy has been omitted. In Hartree-Fock 
calculations |l3| ] these quantum corrections to the ground 
state energy of neutron matter are obviously automati- 
cally incorporated. The Hartree-Fock calculations per- 
formed so far were limited to "spherical Wigner-Seitz 



cells" , which is arguably a reasonable approximation for 
the "nuclei in neutron gas" phase only. To our knowledge 
there exist only one study on this subject where the shell 
effects due to the bound nucleons only however (mainly 
protons) have been taken into account |l4|]. It was de- 
termined that the shell correction energy is smaller than 
the energy difference between different phases and it was 
thus concluded that quantum corrections to the ground 
state energy will not lead to any qualitative changes in 
the sequence of the nuclear shape transitions in the neu- 
tron star crust. 

Our goal is to reach a comprehensive understanding of 
the so called shell correction or Casimir energy in neutron 
stars. There is no well established terminology for the 
energy corrections we are considering here, even though 
the problem has been addressed before to some extent by 
other authors. In the case of finite systems, the energy 
difference between the true binding energy and the liquid 
drop energy of a given system is typically refered to as 
shell correction energy. In field theory a somewhat sim- 
ilar energy appears, due to various fluctuation induced 
effects and it is generically referred to as the Casimir 
energy pf: 
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where go(s) is the density of states per unit volume for 
the fields in the absence of any objects, g(e, 1) is the den- 
sity of states per unit volume in the presence of some 
"foreign" objects, such as plates, spheres, etc., and 1 is 
an ensemble of geometrical parameters describing these 
objects and their relative geometrical arrangement. A 
similar formula can be written for neutron matter energy 
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with the notable difference in the upper integration limit. 
In the above equation go(s, 1) stands for the Thomas- 
Fermi or liquid drop density of states of the inhomoge- 
neous phase and g(e, 1) for the true quantum density of 
states in the presence of inhomogeneities. The parame- 
ters: fi and no are determined from the requirement that 
the system has a given average density 
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Since in infinite matter the presence of various inhomo- 
geneities does not lead to the formation of discrete levels, 
one might expect to refer to corresponding energy correc- 
tion for neutron matter as the Casimir energy. In Ref. 
|jl4| the authors computed a somewhat different quan- 
tity however, than the one we are interested in this work, 
the correction to the ground state energy arising from 
existence of almost discrete levels inside a nucleus in an 
infinite medium. Strictly speaking these levels are not 
discrete, but form narrow energy bands due to the tun- 
neling between neighboring nuclei. The effects we shall 
consider here arise from the "outside" states, which is in 
complete analogy with the procedure for computing the 
Casimir energy. As we shall show, these energy correc- 
tions, arising from the existence of these truly delocalized 
states, are larger than those considered in Ref. iQ. We 
have considered similar issues earlier in finite systems and 
to some extent in infinite 2-dimensional systems as well 
in Refs. 

In order to better appreciate the nature of the prob- 
lem we are addressing in this work, let us consider the 
following situation. Let us imagine that two spherical 
identical bubbles have been formed in an otherwise ho- 
mogeneous neutron matter. For the sake of simplicity, we 
shall assume that the bubbles are completely hollow. We 
shall sidestep the question of stability of each bubble for 
the moment and assume that they are stable and rigid as 
well. We shall ignore the role of long range forces, namely 
the Coulomb interaction in the case of neutron stars, as 
their main contribution is to the smooth, liquid drop or 
Thomas-Fermi part of the total energy. Under such cir- 
cumstances one can ask the following apparently innocu- 
ous question: "What determines the most energetically 
favorable arrangement of the two bubbles?" According to 
a liquid drop model approach (completely neglecting for 
the moment the possible stabilizing role of the Coulomb 
forces) the energy of the system should be insensitive 
to the relative positioning of the two bubbles. A similar 
question was raised in condensed matter studies, concern- 
ing the interaction between two impurities in an electron 
gas. In the case of two "weak" and point-like impurities 
the dependence of the energy of the system as a function 
of the relative distance between the two impurities a is 
given by (spin coordinates are not displayed) 



where kp is the Fermi wave vector and m is the fermion 
mass 
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This asymptotic behavior is valid only for point-like im- 
purities, when kpR "C 1, and where R stands for the 
radius of the two impurities. This condition is typically 
violated for either nuclei embedded in a neutron gas or 
bubbles, when typically kp R 1. As we shall show, in 
the case of large "impurities" (when kpR ^> 1) the inter- 
action energy changes in a rather dramatic manner. If 
one replaces the "weak" impurities with "strong" point- 
like impurities, only the magnitude of the interaction 
changes at large distances, but not the form |17|. The 
interaction (^|) has a pure quantum character, and any 
"noise" (e.g. temperature) leads to a quick disappear- 
ance of the oscillatory behavior and with it of the power 
law character, and the regular Debye screening (which is 
exponential in character) sets in instead. 

The lesson one can learn from this analysis however, 
is that quantum corrections are most likely responsible 
for the interaction of two bubbles/nuclei embedded in a 
Fermi gas and the form of the interaction (|J) suggests the 
most natural way to proceed. The argument of the cosine 
is nothing else but the classical action in units of fi of the 
bouncing periodic orbit between the two impurities. The 
exact form and magnitude of the coefficient in front of 
the cosine can be obtained in a semiclassical approxima- 
tion only after a careful estimation of the leading order 
correction to the leading semiclassical result. Using the 
3-dimensional extension of the semiclassical approxima- 
tion to the so called small disks problem p0| ], we were 
able to obtain a significantly simpler and more transpar- 
ent derivation of this interaction than the original deriva- 
tion JljJ as follows. The correction to the single-particle 
propagator, which depends on the presence of the two 
weak widely separated point-like impurities (R/a <C 1 
and kpR <C 1) is 

5G(r,r',k) cx G (r, n, fc)G (ri, r 2 , fc)G (r 2 , r', k) 

+ G (r, r 2 , fc)G (r 2 , n, fc)G (r!, r', k), (7) 



where 
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where x( r l — r 2 — a ) is the Lindhard response function 
of a homogeneous Fermi gas and Vi(ri) and U 2 (r 2 ) are 
the potentials describing the interaction between impuri- 
ties and the surrounding electron gas. At large distances 
kpa 3> 1 this expression leads to the interaction first 
derived by Ruderman and Kittel f|l8|Jil 



is the free single-particle propagator. Since only "peri- 
odic orbits" contribute to the density of states, the cor- 
rection to the density of states, due to the presence of 
the two impurities and which depends on their relative 
separation only is given by 
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and the corresponding correction to the ground state en- 
ergy is given by the obvious formula 

5£(\ti -r 2 |) oc J kdk - fij Sg(k,\r 1 -r 2 \) 

oc (1Q) 
|ri - r 2 | d 

Only the leading term in the limit |ri — r 2 | — ► oo is 
explicitly shown here. The proportionality coefficient is 
naturally determined by the impurity strength. 

The formation of various inhomogeneities in an other- 
wise uniform Fermi gas can be characterized by several 
natural dimensionless parameters, kpa 3> 1, where as 
above a is a characteristic separation distance between 
two such inhomogeneities, kpR 3> 1, where R is a char- 
acteristic size of such an inhomogeneity, and kps w 1, 
where s is a typical "skin" thickness of such objects. The 
fact that the first two parameters, kpa and kpR, are both 
very large makes a semiclassical approach natural. Since 
the third parameter, kps, is never too large or too small, 
one might be tempted to discard a semiclassical treat- 
ment of the entire problem altogether. However, there is 
a large body of evidence pointing towards the fact that 
even though this parameter in real systems is of order 
unity, the approximation kps <C 1, which we shall adopt 
in this work, is surprisingly accurate Moreover the 
corrections arising from considering kps = 0(1) should 
lead to an overall energy shift mainly, which is largely 
independent of the separation among various objects em- 
bedded in a Fermi gas. On one hand, this type of shift 
can be accounted for in principle in a correctly imple- 
mented liquid drop model or Thomas-Fermi approxima- 
tion. On the other hand, the semiclassical corrections 
to the ground state energy arising from the relative ar- 
rangement of various inhomogeneities have to be com- 
puted separately, as they have a different physical na- 
ture. We are thus lead to the natural assumption that 
a simple hard-wall potential model for various types of 
inhomogeneities appearing in a neutron Fermi gas is a 
reasonable starting point to estimate quantum correc- 
tions to the ground state energy, see Refs. pr^ , p7| , ^l| , p2[ 
and earlier references therein. We shall refer to these 
quantum corrections to the energy as shell effects in the 
rest of the paper. One might expect that such simplifi- 
cations will result in an overestimation of the magnitude 
of shell effects, but the qualitative pattern should remain 
the same. 

We shall consider spherical bubble-like, rod-like and 
plate-like phases only here and we shall estimate the shell 
correction or Casimir energy arising due to a regular ar- 
rangement of such inhomogeneities in an otherwise ho- 
mogeneous neutron gas. One can distinguish two types 
of "bubbles": i) nuclei-like structures embedded in a 
neutron gas and ii) void-like structures. By voids we 
mean the regions in which the nuclear density is signifi- 



cantly lower than in the surrounding space. In the first 
case i), the single particle wave functions can be sep- 
arated into roughly two classes, those localized mostly 
inside the nuclei-like structures and those which are com- 
pletely delocalized. A fermion in a delocalized state will 
spend some time inside the "nuclei" too, but since the 
potential experienced by a nucleon is deeper there, the 
local momentum is larger and thus the relative time and 
relative probability to find a nucleon in this region is 
smaller. One can approximately replace then the "nu- 
clei" with an effective repulsive potential of roughly the 
same shape. In the case of a "bubble" , when the prob- 
ability to find a nucleon inside a "bubble" is reduced, 
again such an approximation appears as reasonable. The 
"nuclei" and "bubbles" we are refering to here, are not 
necessarily spherical, but could have the shape of a rod 
or plate as well. There are of course a number of "res- 
onant" delocalized states, whose amplitude behaves in a 
manner just opposite to the one we have described here. 
However, the number of such "resonant" states is typi- 
cally small and we thus do not expect large effects due 
to them. Moreover, since such states are concentrated 
mostly inside a "nucleus" or a "bubble" one does not 
expect them to affect in a major way the relative po- 
sitioning of two "nuclei" or two "bubbles" . Neverthe- 
less, these are some issues, which certainly deserve more 
scrutiny in the future, even though we hardly expect that 
a more comprehensive analysis will lead to qualitative 
changes of our conclusions. In all these phases the shell 
effects depend on the structure and stability of periodic 
orbits in the system pll . Except for the plate-like nuclei 
phase, where the shell energy can be computed exactly, 
for other geometries one should calculate the contribu- 
tion from all periodic orbits. This is rather tedious task, 
since they proliferate exponentially as a function of their 
length p3| , and moreover this is not really necessary to 
perform. If one is interested in the gross structure only 
of the shell effects, the contribution of the shortest pe- 
riodic orbits should suffice for defining the gross shell 
structure. (We remind the reader, in an infinite medium 
there are really no shells as in a finite system, but we re- 
fer to the corresponding effects in this manner only, due 
to their similar origin because of the appearance of peri- 
odic orbits.) Since the contribution of any given periodic 
orbit leads to an oscillatory contribution to the density 
of states any suitably chosen energy averaging over the 
spectrum, and in particular a finite temperature as well, 
will leave only the contributions due to the shortest peri- 
odic orbits. Since a periodic orbit of length L will lead to 
detail on an energy scale of the order AE — h 2 n 2 /2mC 2 , 
performing an averaging over an energy interval AE will 
effectively mask the contribution of orbits of length £ or 
larger. Moreover, since the geometry of the rod-like and 
spherical phases admit only unstable (hyperbolic) orbits, 
the longer the orbits, the lesser their contribution is, due 
their decreased stability. 
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The simplest system consist of plate-like nuclei 
the neutron gas filling the space between slabs. The 
energy for this system per unit volume can be easily 
uated: 
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where the exact and the Weyl (smooth) energy 
per unit volume are given by 
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In the above formula 

N = Int 



kpa 



(13) 



(14) 



stands for the integer part of the argument in the square 
brackets, and a = L — 2R is the distance between slabs 
and R is the half of the width of the slab. Here L 3 is the 
volume of an elementary (cubic) cell and the factor '2' in 
front stands for the two spin states. The average matter 
density (the number of neutrons per unit volume) po and 
the smoothed density pweyi are determined by relations 
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Using these formulas one can show that the shell correc- 
tion energy has the behavior 
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where G(x) is an approximate periodic function of its 
argument, for x > 1), G(x + 1) w G(x), with properties 
G(x = n/2) ss and approximately —1 < G(x) < 1. 
Furthermore 
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is the actual density of the neutron gas between the two 
slabs and 



1 



L-2R 
L 
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is the filling factor, which is the ratio of the occupied 
volume to the volume of the cell. One can show also that 
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where pweyi is the Weyl approximation to the density 
and F(x + 1) ~ F(x) is an approximate periodic func- 
tion of its argument too, for x > 1, with properties 
— 1 < F(x) < 0.5 and F(x — n) = — 1. This periodicity 
leads to the clear pattern of "valleys" (fcj?a = (n + 3/4)7r) 
and "ridges" (kpa w (n+l/4)7r) in the profile of the shell 
energy shown in Fig. la. These features of the energy 
and density are naturally related to fact that these quan- 
tities are almost periodic functions in the classical action 
along the only periodic orbit in the system, i.e. in the 
variable S — 2kpa. 

In the case of rod-like and spherical voids we shall use 
the semiclassical theory in order to compute the shell 
energy. Since we are interested only in the "gross shell 
structure" we have to take into account a few of the short- 
est periodic orbits among the nearest neighbors only. The 
lengths of the shortest periodic orbits depend on the lat- 
tice type. In the following we will assume the simple cu- 
bic and simple square lattices for spherical and rod-like 
phases respectively. The expression for the shell energy 
density and the neutron density reads: 
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de, (21) 



where g s heii(£,Li) denotes the contribution to the level 
density due to the orbit Li and gweyi is the smooth level 
density determined using the Weyl prescription Q . 

For the rod-like phase we took into account four orbits 
of the length 2Li = 2(L — 2R) and four orbits of the 
length 2L 2 = 2(L^/2 — 2R). Introducing longer orbits did 
not lead to noticeable changes in the patterns presented 
here. Hence the shell energy per volume is equal to: 
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where A\ = A2 = 4 and the chemical potential /1 is 
determined by the condition: 
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A periodic orbit of the type considered by us gives actu- 
ally a contribution with a factor 1/2, since only half of 
it belongs to a particular elementary cell. Because there 
are two spin states, and thus eight orbits in total, each 
type of orbit eventually is weighted by four. The density 
of states was evaluated using the convolution of the exact 
1 -dimensional density of states and the density of states 
given by Gutzwillcr trace formula for the 2-dimcnsional 
system of disks, which is the cross section of the rod- 
like system we are interested in. In some cases such a 
procedure can lead to spurious contributions, which are 
however rather easy to single out, see Refs. p5| . For a 
given periodic orbit of length 2L^, the shell correction to 
the density of states is given by the following expression: 
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where the summation is over repetitions of this orbit and 

h 2 k 2 , , 

e = -T. — • (25) 
2m v ; 

When one is interested in the gross shell structure then 
the contribution of long orbits as well as the contributions 
due to repetitions of short primitive orbits vanish under 
energy averaging. 

The explicit form of the shell energy and of the fluctu- 
ating part of the density reads: 
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The parameter Kj determines the stability of the orbit 
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The shell energy as a function of the anti-filling factor 

ttR 2 

(relative void volume) u — and po is shown in Fig 

lb. The shell energy has a smaller amplitude then in the 
case of the plate-like phase. This comes about because 
the periodic orbits are now hyperbolic in the plane per- 
pendicular to the rods. Note however that the pattern of 
"valleys" and "ridges" looks very similar to the one for 
the slabs. This is because the main contribution due to 
the classical orbit of length 2Li is the same. There are 
small interference effects caused by the orbit of length 
2L2 however. Since it is longer, this second trajectory 
contributes with a smaller weight. 

For the case of spherical voids there are 26 peri- 
odic orbits between nearest neighbors of three differ- 
ent lengths 2Li = 2{L - 2R), 2L 2 = 2(L\/2 - R) and 



2L 3 = 2(W3- R)- Thus the shell energy and density 
are equal to: 
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The contribution due to one periodic orbit to the fluctu- 
ating part of the level density reads: 
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(24) Hence we get: 
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where A\ = 6, A2 — 12, A3 — 8 respectively. The shell 
energy for the spherical phase is shown in Fig lc. In 

AttR 3 

this case the anti-filling factor is given by u = —. 

3 L 3 

A stronger interference pattern due to the orbits L2 and 
L3 can be seen. The amplitude of the shell effects is 
also lower due to the greater instability of the orbit on 
one hand, and due to the smaller relative volume of the 
scatterers on the other hand. 

In a similar manner one can obtain the interaction en- 
ergy between two isolated bubbles at large separations 
(a = L - 2R > R) 



Eqo 



h 2 kpR 2 cos(2fci7d) 
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When compared with the interaction (J5J) one observes a 
similar behaviour, even though now the two "impurities" 
are large kpR 3> 1. It can be shown however that if 
one computes instead the same energy for fixed chem- 
ical potential, instead of particle number as was done 
here, the bubble-bubble interaction will decay inversely 
proportional to the square of the separation p6| . In a 
recent paper ]27j the Casimir energy for similar arrange- 
ments has been calculated using the semiclassical approx- 
imation. In the case of Casimir energy the situation is 
somewhat simple, since instead of two independent di- 
mensionless parameters, k F R and k F L, only one dimen- 
sionless parameter exists, R/ L. Thus the Casimir energy 
for two spheres has naturally the form 
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with an unknown function F(x). A similar, but much 
stronger result can be obtained for the critical Casimir 
energy p8| , where one can show that the theory is con- 
formal invariant. The authors of Refs. |2^| provide also a 
very compelling argument why the semiclassical approxi- 
mation should be particularly accurate for the calculation 
of the Casimir energy in case of ideal metallic boundaries 
and they show that using only the single periodic orbit 
the Casimir energy for two spheres is given by 
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and dismiss this result as being valid for large sepa- 
rations, since it contradicts their expectations that it 
should agree with the Casimir-Polder interaction [E9l 
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The authors of Ref. j27j argue that the contributions 
arising from the diffractive paths discussed in Refs. [j30[ , 
should eventually lead to additional contributions, which 
will cancel exactly this longer range interaction and in the 
end, the authors hope that the Casimir-Polder result will 
be retrieved. The difference between these two results for 
the Casimir energy is very similar to the difference be- 
tween the interaction (|^) between two point-like impu- 
rities (kpR <§C 1) and the interaction between two "fat" 
bubbles [kpR ^> 1) ^6|. In the case of "fat" bubbles, the 
contribution of diffractive orbits are exponentially small 
(oc exp(— akpR), where a is of order unity) |3(i| ], as one 
would naturally expect in the case when rays are a very 
good approximation to the wave phenomena. The reso- 
lution of this apparent conundrum lies in resolving the 
clear clash of limits. When the size of the scatterer R 
decreases the contribution of the diffractive paths (creep- 
ing orbits) increases and an increasingly larger number of 
them contribute significantly to the scattering and thus 
to the propagator. In the limit kpR — > the standard 
geometric orbit approach has to be modified, see Ref. 
p0| and our discussion around Eqs. (7-10). It is notable 
that in the case of the critical Casimir effect, even longer 
range interactions (oc l/a 1+e , with very small e) between 
two spheres are possible p8| . 

The structure of the shell energies shown in the Fig. 
1 indicates the existence of the optimal void sizes (with 
respect to the shell effects) for a given outside nucleon 
density. Note that for all phases and for po > 0.05/m -3 
the shell energy exhibits a remarkable softness toward 
adding additional neutrons to the system (the "valleys" 
and "ridges" are almost horizontal in the Fig. 1). Hence 
one can conclude that once the size of the voids have been 
determined by minimization of the total energy of the 
system, an increase in the number of neutrons outside 



the voids will not affect much the shell energy of the 
system. However, the surface energy will be affected. 

In the Fig. 2 we show the shell energies as a function 
of po I0r the optimal filling factors and nuclear radii de- 
termined in Ref. || . One can see that the amplitudes of 
the shell energies in the region po ps 0.04 — 0.07 are of the 
order of 10keV/ fm 3 , ikeV/fm 3 and O.ObkeV/ fm 3 for 
plate-like, rod-like and bubble-like phases, respectively. 
There are usually one or two shallow shell energy min- 
ima for the density range po > 0.03/m -3 . The minima 
are more pronounced in the case of spherical bubble- 
like phase mainly due to the stronger interference effects 
caused by longer orbits. 

Once a phase is formed there is a positional order main- 
tained by the Coulomb repulsion between spherical nu- 
clei, rods or slabs |^|-|],|]|7||j| . Although the Coulomb en- 
ergy is a smooth function of the void displacement pi] , 
the shell energy is not. Since several different orbits con- 
tribute to the shell effects (except for slab-like phase) 
the displacement of a single bubble-like or rod-like void 
from its equilibrium position in the lattice will give rise 
to the interference effects. The interference pattern will 
depend on the type lattice. For the simple cubic and sim- 
ple square lattices for spherical nuclei and rods, respec- 
tively, we show in the Fig. 3 the changes in the energy 
due to such "defects" . For the plate-like system there is 
only one direction of displacement (we do not consider 
the shear mode) denoted by x perpendicular to the slab 
(Fig. 3a). Since the rod-like phase is a two-dimensional 
system, in the Fig 3b we have shown the shell energy as 
a function of two perpendicular displacements x and y. 
They are perpendicular to the rods and point in the di- 
rection of the nearest neighbor. The same axes have been 
chosen for the spherical system although it will not ex- 
haust all possible directions in the system. The behavior 
of the shell energy in this case is shown in Fig 3c. 

The structure of the shell energy surface as a function 
of a displacement depends on the lengths of the shortest 
periodic orbits. Except for the trivial plate-like phase, 
in the rod-like and spherical phase there exist directions 
into which is easier to locally deform the lattice. 

In Fig. 4 we show the pattern of the energy changes 
induced by deforming the rod-like lattice. We consid- 
ered only volume conserving deformations. The square 
lattice was stretched by a factor a in the cc-direction, by 
a factor (3 in the y-direction and also the angle between 
the two axes has been changed to 7. In order to preserve 
the volume all these three parameters should satisfy the 
condition 

a/3sin7 = l. (38) 

The case a = @ and 7 = 7r/3 correspond to a perfect 
triangular lattice. 

Increasing the temperature will weaken the shell ef- 
fects. At sufficiently high temperatures the nuclear lat- 
tice will disappear. At smaller temperatures however, 



G 



when the lattice can be regarded as frozen, the rise of 
the temperature will affect mainly shell effects in the neu- 
tron gas. In order to wash out completely the shell effects 
the temperature T should be of the order of half of the 
distance between shells The spacing between two consec- 
utive shells is determined by the length of the shortest 
orbit, a = L — 2R. Thus the energy distance between 
shells can be determined from the requirement: 2ka = 2n 
and is given by the expression: 



2 ^2 



AE 



2m (L -2R) 2 ' 



(39) 



For the optimal filling factors and lattice constants of var- 
ious phases obtained in Ref. El one obtains the following 
estimates for the critical temperature: 



T c w 32MeV for plate-like system, 
T c w 19MeV for rod-like system, 
T c sa 12MeV for spherical system. 



(40) 



An accurate description of the shell effects as a function 
of the temperature can be obtained using the tempera- 
ture averaged level density 0,pl| : 



shell 



^ sinh r iin (T) 



9sheu(£,Li), (41) 



where the sum is taken over all periodic orbits including 
the number of repetitions n of the orbit and 



2nTmnLi 
Wk 



(42) 



Consequently the oscillating part of the free energy 
density is given by the formula: 



F, 



hell 



L 3 
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L3 



(e - fi)g S heii (e,T)de. 



(43) 



The estimates for different phases are shown in Fig. 5. 
For simplicity we have retained in these calculations the 
value of the Fermi momentum kp equal to its zero tem- 
perature limit, and therefore the neutron matter density 
is not temperature independent in these figures. One 
can see that thermal effects will wash out the shell cor- 
rection energy at temperatures of the order of 10 MeV 
and higher. 

Now at the end of this analysis we suspect that there 
are a lot of other effects, which might be relevant. We 
did not consider periodic orbits bouncing between three 
or more objects. An orbit bouncing between two bod- 
ies leads to a pairwise interaction. Orbits bouncing be- 
tween three or more bodies would lead to genuine many 
body interactions. We have also considered only per- 
fectly smooth objects. If one allows for some degree of 
corrugation of these surfaces, many more periodic orbits 
are likely to appear and that would lead to even more 



complicated interactions and more complicated interfer- 
ence patterns. The fact that corrugation can influence 
in a significant, perhaps major way, the Casimir energy, 
has already been predicted and measured experimentally 
p2f . The long range character of the interaction together 
with its oscillatory nature could very easily be at the 
origin of disorder, even at zero temperature. At finite 
temperature disorder is more likely to occur, due to en- 
tropic effects JTfjJlTj • We did not consider here the role of 
pairing, which we expect however to lead to a certain flat- 
tening of the shell effects [jL6|]l7|] , which, however, should 
not be interpreted as disappearance of shell effects. Espe- 
cially at subnuclear densities neutron pairing should be 
rather strong [B3[ . A completely different type of softness 



of these structures has been argued in Ref. 34 1, accord- 
ing to which the mantles of neutron stars resembles more 
liquid crystals than solids. 

In the paper we have studied the shell effects in the 
neutron medium filled by different nuclear phases. To 
our knowledge this is the first approach which considers 
specifically the shell effects in the outside neutron gas and 
we aimed at discussing its basic features. Even though 
in principle Hartree-Fock calculations include in princi- 
ple such effects already, the calculations performed so far 
|]l3f were too narrow in scope and did not address this 
issue specifically. Using semiclassical methods, we have 
analyzed the structure of the shell energy as a function of 
the density, filling factor, lattice distortions and temper- 
ature. We expect that our result overestimate somewhat 
the amplitude of the shell effects. However, the emerging 
qualitative overall picture should remain valid and fur- 
ther microscopic studies are highly desirable. The main 
lesson one should remember from this work is that the 
amplitude of the shell energy effects is comparable with 
the energy differences between various phases determined 
in simpler liquid drop type models. The magnitude of 
the quantum corrections to the ground state energy of 
the inhomogeneous neutron matter we have found is sig- 
nificantly larger than that determined in Ref. |l4j]. The 
analysis of Ref. [O] was limited however to the motion of 
nucleons inside nuclei embedded in a lower density neu- 
tron gas. 

Our results suggest that the inhomogeneous phase 
has perhaps an extremely complicated structure, maybe 
even completely disordered, with several types of shapes 
present at the same time. 
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